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Abstract. Definition of generalized normal form for a system of ODEs corresponding 
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1. Introduction 

In the introduction, some necessary facts from symplectic and contact geometry are 
given. For a complete overview see [TJ Ch 1, §§ 6, 7] and references therein. 

A symplectic form on a manifold M of even dimension s = 2 n is a closed 2- form 
u of maximal rank, i. e. du = and u n ^ everywhere on M. There exists a local 
coordinate system zi, . . . , z 2n for which u = ^2i =l dz n+i A dzi (see [U p 154]). Such 
coordinates are called canonical coordinates. 

A transformation / (or an infinitesimal transformation X, respectively) of M is 
said to be symplectic if f*u — u (L^u = 0), where /* denotes the linear mapping 
transposed to df, and Lx is the Lie derivative along the vector field X. 

In canonical coordinates, an infinitesimal symplectic transformation is given by 

where H is an arbitrary differentiable function. H is called the Hamiltonian, and 
infinitesimal transformation ([I]) is called a Hamiltonian vector field. The corresponding 
system of ODEs is called a Hamiltonian system. 

Defining the Poisson bracket of two functions $ and \I/ by Xi^ t m = [X$, X^], in 
canonical coordinates we have the classical formula 

{$, = \ — 

A contact form on a manifold M of odd dimension s = 2 n + 1 is an open cover 
{U a } of M together with a system of 1-forms a; Q of maximal rank, i. e. u a A (du a ) n ^ 
everywhere on U a , such that u a = (p aj 3 up on U a C\Up for some functions (p a p ^ 0. 

For any point p G U a , there exists a coordinate neighbourhood U p C C/ a with 
coordinates z\, . . . , ^n+i in which u a = dz 2n+ i — Yji=i z n+idzi (see [TJ p 150]). Just as 
in the symplectic case, we call such coordinates canonical. 

A transformation / of M is said to be contact if for all a, (3, f* u a = A a p up, where 
A a/ 3 7^ everywhere on f~ l (U a ) C\Up. Accordingly, an infinitesimal transformation X is 
called contact if Lx u a = X a u a for some functions A a . 

In canonical coordinates, an infinitesimal contact transformation is given by 



x - y I" dH d i ( dH i z dH ) 9 

^ dz n+i dzi \dzi n+t dz 2n +iJ dz n+i j 



+ 



where if is an arbitrary differentiable function. The function H is sometimes called a 
contact Hamiltonian, or simply the Hamiltonian. 
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In the contact case, the Poisson bracket in canonical coordinates has the view 

{$, ^} = $- — + 

OZ2n+l OZ2n+l 

£~f\dzi n+l dz 2n +i J dz n+i \dzi n+l dz 2n +i ' dz n+i ' 

From now on, only functions that are analytic in some neighbourhood U of the 
origin in R s are considered. We denote the algebra of analytic functions in U by WLu(z), 
and the algebra of polynomial functions in z%, . . . , z s by M.[z]. 

Diffeomorphisms of M s onto itself of the form exp(Xp) with P 6 M.[z] generate 
the subgroup 2l s in the group of symplectic (if s = 2n), or contact (if s = 
2n + 1), transformations of M. s . The algebra of infinitesimal symplectic, or contact, 
transformations is closed under the action of this subgroup, and the Hamiltonian 
functions are transformed according to the law 

H ->■ exp(P)(if), 

where P( ■ ) = {P, ■ }. 

Let the origin be a singular point of the vector field X^, i- e. Xh\ z=0 = 0, and let 
the corresponding Hamiltonian have the view 

H = H + H l (H G R[z], H x e Ru(z)). (4) 

We call H the unperturbed part, and Hi the perturbation of the Hamiltonian H. 

In the first four sections of the present article, the general problem of bringing the 
perturbation of H into the simplest form (normal form) by transformations from the 
group 2t s that do not change the unperturbed part of H is being solved. Such problems 
often arise in studying Hamiltonian and nearly Hamiltonian systems. As a result, there 
are many different ad hoc methods and definitions of normal forms (see, for instance, 
[21 El H]) that require the unperturbed Hamiltonian to satisfy certain conditions such 
as quadraticity, homogeneity, etc. But such requirements can be excessive in some 
physical problems. The generalized normal form introduced in this article is free of 
such restrictions. It is worth noting that the contact case was earlier considered only in 
the simplest, generic, situation. However, it turns out that real physical problems, 
for instance thermodynamics of plasmas, often refer to the non-generic cases. In 
section 5 we introduce some applications of the method of generalized normal forms 
to thermodynamics of non-ideal media. 



2. Basic definitions 



As usual, the sets of positive and nonnegative integers will be denoted by N and Z + 
respectively. 

Let neN and s = 2 n, or 2 n + 1, be fixed. 
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Definition 1. A vector 7 = (71, . . . , 7 S ) G N s is called a weight of the variable 
z = (z 1 , . . . , z s ) if GCD (ji, . . . , 7 S ) = 1. If, moreover, for some natural a > 2, 
71 + 7„ + i = . . . = 7„ + 7 2n = <7, and in the case of odd s, 72«+i = cr, we say that 
7 is a canonical weight. 

Definition 2. Given a power series (or a polynomial) P = XLez s a ^ the generalized 
order (the generalized degree) of P with weight 7 is the least (respectively, the greatest) 
of numbers v ■ 7 = 1/1 71 + . . . + v s 7 S such that a v 7^ if P ^ 0, and 00 (or 0, respectively) 
otherwise. 

Definition 3. A polynomial P G M[z] is said to be quasi-homogeneous with weight 7 
if either its generalized order and generalized degree with weight 7 are equal or P is 
identically zero. 

Definition 4. Given an arbitrary power series P = E„ eZ ^ a vZ v , define its [&]-jet with 
weight 7 by jif ] (P) = E* 7 =o a ^- 

Denote the generalized order of P with weight 7 by ord 7 P, the generalized degree 
by deg 7 P, and the vector space of all [A;] -jets with weight 7 by 0^- 

Consider a power series H(z) = Hq{z)+H\{z) with the unperturbed part Hq G 
and the perturbation Px G such that for some fixed canonical weight 7, 

ord 7 Hi > deg 7 Po > ord 7 Po- (5) 

For brevity, denote 5 = ord 7 Hi — ord 7 P and A = ord 7 Hi. 
Define an inner product on polynomials by 

((P, Q» = P(9) Q(z)\ z=0 (P, Q G R[z], d = (d/dzi, . . . , d/cb s )), (6) 

and fix some linear order >- on the set of monomials in z. 

For any integer d > A, define a vector space 91^ C 1 ' as follows: 

0t rf = {jif +<5 ~ 1] ({P, P }) : FeRN, ord 7 P>rf + 5- A + a}, (7) 
and denote the set of leading monomials of the elements of 9T rf with respect to y by 

lm (m d ). 

Lemma 1. For any integer d > A, there exists a basis <& d = {Gi} of 91^ such that 
{{Gi, LM (Gj))) = 5ij, where 8^ denotes the Kronecker delta. Herewith, the sets LM (0 d ) 
and LM (9Lj) are i/ie same. 

Proof. The proof is by induction. Let {Pi} be an arbitrary basis for 91^. Setting 
d = Pi/((Pi, LM (Pi))) gives LM (Gi) = LM (P x ) and ((d, LM (d))) = 1. 

Suppose that for some k G N and for all i, j < k, we have ((Pj, LM (Pj))) = 
Denote G" fe+1 = P fc+1 - £ti« s *+i> LM ( s 0)> s i> so that 1 G 'k+n LM (P,))) = 0. Then 
setting G k+ i = G' k+1 /{{G' k+v LM (G' k+1 ))) and G t = B t - ((P 4 , LM (G fc+1 )» G fc+1 for all 
i < k gives LM (G-) = LM (P;), and hence {{G h LM (G,))) = 5, tj for all i, j < k + 1, 
which completes the induction and proves the first statement of the lemma. 

Since 9Lj is finite-dimensional, the second statement follows at once from the fact 
that all of the elements of &d have different leading monomials. □ 
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Definition 5. Define a minimal resonant set of index d as a complement of LM (Ot^) 
in the set of monomials of generalized degree k, where d < k < d + 5. 

Denote the vector space spanned by a minimal resonant set of index d by Then 
given a basis &d for 0^ as in Lemma 1, the formula 

n d (P) = P- Y,((P, LM (Gt))) d (8) 

i 

defines a projection of the vector space of [d + 5 — l]-jets of generalized order greater 
than or equal to d onto 

3. Normal form theorem 

Theorem 1. For any integer N > A, there exists a transformation from 2l s that brings 
H into the form 

l(N-A)/6] 

H = H + R k + S > (9) 

k=0 

where Rk G ^A+kS and ord 7 S > N. 

Proof. Let F G M[z] and ord 7 F > 5 + o~. Then the transformation exp(X^) G 2t s brings 
H into the form 

exp(F)(H) = H + H 1 + {F,H } + ..., 

where ord 7 ({F, H }) > A, and the dots denote terms of generalized order greater than 
or equal to A + 5. 

More generally, if ord 7 F > k 5 + o for some k G N, then in the same formula 
ord 7 ({F, H }) > A + {k — 1) 5, and the generalized order of terms denoted by dots is 
greater than or equal to A + k 5. 

Suppose that H = H + ^^Zq 1 Rk + S m , where Rk G 9^A+fc<5, ord 7 S m > A + m5, 
for some m G Z + . The proof is by induction on m with So = H\ for m = 0. 

Write S m as S m = P m + Q m , where P m = J 7 A+( - m+1 ^ 5 ^(5 m ). Since P m — 
^A+ms{Pm) G yiA+mS, where 7r A+m 5 and 9Ta -\-m8 defined in (GO)-®, there exists 

a polynomial F m of generalized order greater than or equal to (m + 1) 5 + a such 
that P m + J [ f +{m+1)S ~ 1] ({F m} Ho}) = Tt A+m s(P m ). Denote R m = 7r A+m5 (P m ). Then the 
transformation exp(Xi? m ) G 2l s brings if into the form H + YUk^^ + Sm+i, where 
ord 7 S m+ i > A + (m + 1) 5. 

By induction, m can be made arbitrarily large, and in particular, equal to 
\{N - A)/ 5], so that ord 7 S m+1 > A + ([(N - A)/5] + 1)5 > N. Since at each step 
the normalizing transformation has the form exp(Xp m ) G and the number of steps 
is finite, the composition of these transformations lies in □ 

Definition 6. We say a Hamiltonian H with the unperturbed part Ho, as well as the 
corresponding vector field, to be in generalized normal form up to generalized degree N 
if it has the form (Q. 
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4. The case of quasi-homogeneous unperturbed part 

In case where the unperturbed Hamiltonian is quasi-homogeneous for some canonical 
weight, there is an alternative way of describing generalized normal forms in terms of 
what we call the "resonant sets". This method requires no linear order on monomials 
and is stated as follows. 

Denote the vector space of quasi-homogeneous polynomials in z of generalized 

Ik] 

degree k G Z + with weight 7 over R by RJy . 

Let H = H [x] G where 7 is a canonical weight, and x > °~- Then for each 
k G Z+, the linear operator if 7 x ' = {Hj, ■ } maps K 7 to K 7 x a \ and its conjugate 
with respect to the inner product fl6]) has the form 



i=l 

in the symplectic case, and 
H^f = z 2n+1 H^ - [ dH ^ 



dH [ *\* ^ (8H [x] 

z? 



dzi J V dz n +i 



-dZ2n+l ■ 

^ fdH\ x] 8H [x] y ( 8 \ fdH[ x] 

1 / Z n+ i I — h Z n+ i — I — I Zi + Z2n+1 Ti ) I — 

^ V dzi oz 2 n+i^ v oz n+i J \dz ni 

in the contact case. Here P* = P(d) for each P G M[z], 
Definition 7. We call the equation 

HM*(P) = (10) 
for P G R[z] the resonance equation, and its solutions resonant polynomials. 
Denote the vector space of resonant polynomials by 9 c t 7 . 

By properties of the inner product (EJ), it follows that Diy splits into direct sum of 
the orthogonal subspaces spanned by quasi-homogeneous resonant polynomials of 
generalized degree k G Z+. 

Definition 8. We call a set of quasi-homogeneous polynomials &^ = {SjY^ a 
resonant set of generalized degree m with weight 7 if for some basis {i2j}|™ x of 
det({((i?j, 5^)) }f j =1 ) 7^ 0. Given a resonant set &^ for each generalized degree m G Z 
we call their union, (3 7 = Um=o®7 > a re sonant set. If, moreover, (5 7 consists only of 
monomials, we say that it is a minimal resonant set. 

Clearly, the definition of a resonant set does not depend on the choice of basis 

\k] 

for DVy*. Note that in the case of a quasi- homogeneous unperturbed part, a minimal 
resonant set in the sense of Definition [5] is also a minimal resonant set in the sense of 
Definition [HI and vice versa. 

\y\ 

Theorem 2. For each integer N > A and arbitrarily chosen by H\ resonant set (5 7 , 
there exists a transformation from 2l s that brings H to the form 

N 

H = H + Y,R? + S, 

k=A 



in 

■ 
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where e span (©if') and ord 7 S > N. 

Proof. The theorem follows from Definition [8] and the Fredholm alternative for finite- 
dimensional spaces. Detailed proof in the symplectic case (in the contact case, the proof 
is exactly the same) is given in [5]. □ 

For the symplectic close result to Theorem 2 can be found in [4|. However, 

the normal form equations given there are very complicated, since the unperturbed part 
at each step of normalization involves normalized terms of the perturbation (for further 
information, see [6]). The quasi- homogeneous symplectic case was first analyzed in [5], 
where the canonicity condition on the weight was stated and definition of a resonant set 
was introduced. 



5. Applications of generalized normal forms in thermodynamics 

5.1. Thermal equations of state 

In thermodynamics, relations between macroscopic parameters are described by the 
equations of state. Equations that provide a relationship between the temperature T, 
the volume concentrations rti (i = 1, k) of the system components, where k denotes the 
number of components, the pressure P, and other generalized thermodynamic forces, 
are called thermal equations of state. For example, 

k 

p = tJ2^ (ii) 

i=l 

t=i i=l 

are, respectively, thermal equations of state of the ideal gas and the classical Debye- 
Hiickel plasma written in ESU CGS units, where T is considered in erg (see, for instance, 
[3 pp 151,282]). Here, denotes the electric charge of a particle of i-th. type. 

There are two approaches to describing thermodynamical properties of non-ideal 
media, the rigorous one treating media as a system of nuclei and electrons interacting 
by Coulomb's law, and the semi- empirical one that consists in constructing qualitative 
models, in which the general form of functional relationships is established on theoretical 
arguments, and the numerical coefficients are determined experimentally | s 8], 
Method of generalized normal forms allows to reduce significantly the number of 
unknown coefficients in such qualitative models of non-ideal media. 

Consider a thermal equation of state of the form 

H(P, T, n) = 0, n=(m, ...,n k ). (13) 

Using the Gibbs-Duhem equation (see, for instance, [91 p 322]) 

k 

dP — s dT — ^ dfii = 0, 

i=l 
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where s denotes the entropy density, and //j the chemical potential of the i-th type 
of particles, ffTB"]) is seen to be a first order partial differential equation for P. Its 
characteristic equations are determined by infinitesimal contact transformation (j2J): 

OH 



ni = n, t 
f = 



where subsystem (fT4"1i ) contains only variables that enter (fT3"|) . Here, the dot denotes 
differentiating with respect to some parameter that does not have a physical sense itself. 
Note that the quantities T, nj/^2 i n i , H/J2i n i are integrals of (JHJ), and the surface 
given by ( TT3|) is integral for this system. 

Suppose that H can be represented in the form PJ-©, where Zi = Zk+i = T, 
Zk+i+i = n iy Z2(k+i) = s, Z2(k+i)+i — P (i — 1, k). Then the above theory is applicable, 
so by Theorem 1, H can be put into normal form if up to a desired degree by a contact 
transformation, with the surface H = transformed to H = 0. 

By definition, variables s and \ii do not enter a thermal equation of state. So if we 
do not want them to appear in H, it is necessary to consider the transformations that 
depend only on P, T, and rij. This follows from the form of the Poisson bracket ([31). 

5.2. Perturbation theory for ideal gases 

Equation of state of a mixture of non-ideal gases is usually written in the form of a 
virial expansion, also called the Mayer cluster expansion: 

k 

P = T ^ni+ ^ B rn n m , m = (mi, ...,m fc ). (15) 

i=l |m|>2 

B m are called virial coefficients and characterized by contributions of the interactions 
within |m|-particle groups, or clusters, of molecules of gas. 

Suppose that B m are analytic functions of temperature near the point T = 0. Then 
the following theorem holds. 

Theorem 3. For each integer M > 2, all the virial coefficients B m with \m\ < M in 
expansion (T5}) can be eliminated by a contact transformation. 

Proof. In generalized degree 2 with weight 7 = (1, . . . , 1, 2), where 2 is the weight of 
P, equation of state ( |T5l) can be written as 

i=l ij=l 

Here b is a constant symmetric matrix. Clearly, the contact transformation P = 
P — Yli j=ibij n i n j/2, fii = Hi — Yl k j=i^ij n j brings it into the form (ITTi) . Thus, the 
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part of the second virial coefficients that does not depend on temperature can always 
be eliminated by a contact transformation. 

Let us prove the statement for terms of the perturbation of generalized degree 
greater than 2. For this purpose write out the resonance equation (flOl) for the 
unperturbed Hamiltonian = P — T ^2i=i n i corresponding to (TTTj) . taking into 
account that the perturbation depends only on P, T, and : 

dR A OR 



P dP + ^ ni d^~ R = (-R ^ K[P, T, Hi}). 

i=i 

Monomials of the form P T l and rij T l make up a basis of solutions of this equation, and 
hence by Theorem 2, normal form of the perturbation. In particular, terms of degree 
higher than one in Hi are absent in normal form. □ 

Theorem 3 agrees with a well known fact that the media described by virial 
expansions (for example, the Van der Waals gas) permit continious transitions between 
any two states without ever separating into two different phases [7J pp 298,299]. 
Apparently, for the same reasons as claimed in Theorem 3, the virial expansion turned 
out to be useless for describing critical phenomena, such as the gas- liquid- solid triple 
point (for further discussion, see PHI pp 105-110]). 

5.3. Perturbation theory for classical Debye-Hiickel plasmas 

Equation of state (|T2|) of fc-component plasma holds under assumptions of electrical 
neutrality and smallness of the Coulomb interaction in comparison to the kinetic energy 
pp 279,280]: 

k k _ „ 



i=l i=l 

The unperturbed Hamiltonian corresponding to (fl2|) can be chosen as 

k k 

H = T(P-Tj2n?} 2 -l(j2^ q *) 3 - (17) 

i=l i=l 

Without loss of generality, assume that qi = 0, q 2 > 0, and q 3 < 0. Then the 
change of variables n x = Yli=i n u n 2 = Z)i=i n i% n 3 = (n/9) 1/3 Yl^i^Qh ™j = %> 
where j = 4, k, is nondegenerate and can be extended to a contact transformation. For 
this purpose, choose fii that satisfy d/Ij = rijd/ij {% = 1, k), and leave P, T, and s 
unchanged. Notice that condition ffToTi ) takes the form n 2 = 0. 

Let us denote 

zi = Mi : z S = z* = T ( . = T) (lg) 

Zk+l — ^1, Zk+j — %+l, Z2k — S, Z2k+1 — P 

In new variables, the Hamiltonian (fTTl) is written 



Hq — Zk (z 2 k+i — Zk Zk+i) 2 — zf +2 . 
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Inequality (TTBT? ) . which takes the form z k+2 "C z\ according to (1T51) . implies a 
natural condition on weight 7 of the variable z : 7^+2 = 3 7^. It follows in turn from (fT2l) 
that 72/k+i = 47 fc . Furthermore, ( 1T8|) implies 7^+^ = 7^+1 (j = 2, A; — 1). The remaining 
components of 7 are determined by the canonicity condition. Thus, in the case under 
consideration, the natural canonical weight has the form 7 = (1, . . . , 1, 3, . . . , 3, 4), 
and Hq is a quasi-homogeneous polynomial in z of generalized degree 9 with weight 7. 

Theorem 4. The minimal resonant set & corresponding to a lexicographical order 
determined by the order z 2 k+i y z k y z k+1 y . . . y z 2k -i consists of the following 
monomials: 

ft) Z k ' ' ' Z 2k-1 i 

(ii) z 2k+1 4 fc+1 z k+j, where j = 1, k - 1, and z\ k+l z k k+1 ; 

(Hi) ^ 2 fc+l Z k+1 ' ' ' Z 2k-1 > exce Pt z 2k+l z f +2 , 

where i 2k+1 , i k , ik+i, kk-i € 

Proof. For all i 2k+ i, i k , i k +i, ■ ■ ■ , i2k-i e we have 

J7T *2fc + l ik *2fc-11 _ *2fc+l~l ifc _*2*-l 

l JJ 0) ^2fc+l k ' ' ' ^2Jfc-l J ~~ ^2k+l z k ' ' ' ^2fc-l A 

l -^Y~ + z fc+ i + • • • + iajfe-i - l) z 2 2k+1 z k - (19) 

— 2 (%fl + . . . + »2fe— 1 ~~ 1) ^2fc+l ^fe Zk+l — ^2k+l z \ Z k+1 + 2 i 2k+ i z\+ 2 

This implies that all monomials divisible by z 2k+i z k} except z 2k+i z b k k+1 z k+ j for j = 
1, k — 1 and z 2k+l z k k+1 } do not belong to ©. 

Furthermore, in case where i 2k+ i + 2 z^+i + . . . + 2 «2fc-i = 2, the right-hand side of 
(fl9~|) is identically zero if Z2fc+i = 0. Reducing one by one the leading monomials in the 



Poisson bracket ( !T9|) in the case where i 2k +\ = 2 and i k+ j — (j = 1, k — 1), we come 
to the following equality: 

{#0, (^2fc+l ~ Z k Z k +l) 2 } = 4:Z k k (z 2k+ i Z k+2 — Z k Z k+ i zl +2 ). 

This implies that the monomial z 2k+ i z\ +2 does not belong to (3. For i k > 1, further 
reduction of the leading monomial gives identical zero. 

Collecting all the results together, we come to the statement of the theorem. □ 

From now on, we consider, for simplicity, the plasma consisting of the hydrogen 
atoms, protons and electrons. We denote the corresponding parameteres by the 
subindices a, p, and e. In this case, q a = 0, q p = —q e = e, where e is the elementary 
charge, n e = n p , and 

32/3 

Z\ = // a , z 2 = - — . (/i p + fl e — 2 /i a ), z 3 = T, 

"I 1/3 2 (2°) 

Z4 = n a + 2 rip, z 5 = 32/3 n p , z 6 = s, z 7 = P. 
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Theorem 4 gives the following generalized normal form of the equation of state: 
z 3 (z 7 -z 3 z 4 ) 2 -z 3 + Yl h^^4z{z k 5 + h^°^z\44 + 

i+3j+3fc>10 4i+3j'+3fc>10 

oo 

+ J2[h [3 ' l '°'° ] Z 3 7 + (h^ i+2 ^ Z 4 + /^+2,0,l] Zh) Zl -|] z i = Qj (21) 
i=l 

where ti 1 ' ' ^ = 0. 

Let us consider the case, where the perturbation consists of elements of the minimal 
resonant set & of generalized degree 10. According to (121j) and the physical restriction 
Z7/Z3 z a as either z 3 — > 00 or z^ — y 0, the form of such perturbation is determined 
uniquely up to three constant coefficients 

= - Z5 (az 4 + bz 5 + c zl) (z 7 - z 3 z 4 ). (22) 

Note that if we considered an arbitrary perturbation (not in normal form), we should 
have additionally taken into account the terms z 2 z 2 and z 7 z\. 

In order to find out the physical meaning of a, b, and c, first look at the behaviour 
of the Debye-Huckel plasma (a, b, c = 0) ignoring (fTEk ). By integrating (fll]) . we come 
to the following expressions for z\ and z 2 ■ 



Zl = T In z, + C u z 2 = -3 2 / 3 tt 1 ^ e J + C 2 , 

y (1 + a) T 

where a = n p / (n a + n p ) is called the ionization coefficient. Recall that a is an integral 
of ( 1T4|) . Since at low values of the concentration Z4, plasma turns to a mixture of ideal 
gases, C\ and C 2 can be determined by replacing the chemical potentials of components 
in ( 120|) with the chemical potentials of ideal gases that have the form 

rn /2rf\ 3/2-1 
//id = T In - — j 

where n denotes the concentration, Z the partition function of a particle, and m the 
particle mass p 163]. For electrons and protons Z = 2, and for the hydrogen atoms, 
for simplicity, we take Z = 4 exp(Ry/T), where Ry f» 2.18 • 10~ n erg is the Rydberg 
constant. This partition function ignores all the energy levels of the hydrogen atom but 
the lowest one. For more precision the Planck-Larkin partition function should be used. 
Solve the equation of chemical equilibrium that takes the form 

zi + 3~ 2/3 2 tt 1/3 e 2 z 2 = (23) 

to obtain 



a 



2vr/i 2 m a \3/2 /R y e^STrnp 

e^n 



(24) 



1 + n p [ ) exp 1 — — — 

L V m e m p 1 J \ 1 1*1 A 

which implies the effect of the ionization potential decreasing as n p increases. Note that 
this expression is equivalent to the classical Saha formula with the partition function 
of the hydrogen atom replaced by Z = 4 exp(Ry/T — e 3 ^8 n n p /T 3 ^ 2 ). Such approach 
to perturbation theory of plasmas consisting in perturbing the hydrogen atom partition 
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(a) (b) 

Figure 1: a(n Q ) and P(n ) in the Debye-Hiickel model at T = 20 000 K (— — ), 
T = 23 000 K (- -), and T = 26 000 K ( ). 



function is usually called a chemical model and used, for instance, in [TT] where the 
Planck-Larkin and the nearest neighbour partition functions are used in this context 
(for further information on chemical model, see [T2"l ss 5.1,5.2]). 
Denote no = n a + n v and rewrite (1T21) as follows 

P = Tn (l + a)-^^(an f 2 . (25) 

Then (J25J) and (T25]) together with n p = an give a complete thermodynamical 
description of Debye-Hiickel plasma. 

According to (1241) . in the low density region where (fl6l ) is satisfied, the Saha 
approximation is valid. As the density no increases the ionization coefficient stabilizes 
at first, and then after n exceeds the value ~ 8 • 10 21 particles/cm 3 , rapidly tends to 1 
(see figure [Ta]). At temperatures less than « 23 000 K the curve P(n ) has two branches 
where (dP/dn ) T > (see figure [Tbl) which corresponds to splitting of the plasma into 
two stable phases. The thin one is weakly ionized and called the dielectric gas phase 
and the dense is strongly ionized and called the metal gas phase. This phenomenon is 
called the plasma phase transition. For temperatures greater than « 23 000 K metal 
gas phase becomes unstable since (<9P/<9no)r is negative on the corresponding branch. 
An overview of relevant theoretical results can be found in [T2l s 6.4] and references 
therein. It should be noted that plasma phase transition has not yet been detected 
experimentally. For this reason it is often called hypothetical in literature. 

As can be seen from figure [Tbl the density n is bounded above. Beyond this bound 
the pressure becomes negative, which is meaningless. Apparently, at higher densities 
a thermodynamical description of plasmas is impossible due to absence of quasi-static 
states which are the only considered in thermodynamics. 

System (1141) for perturbed Debye-Hiickel plasmas with the perturbation ( 122]) can 
be integrated explicitly as well. The results for small perturbations are given in figure [2j 
The coefficients a and b affect the pressure and the concentrations of phases at the point 
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Figure 2: a(n ) and P{tiq) in perturbed Debye-Hiickel plasmas at T = 23 000 K for 

(a)-(b): a = ( ), aT/e 4 « -0.04 (- -), and aT/e 4 w 0.04 ( ), 6, c = 0; 

(c)-(d): 6 = ( ), 6T/e 2 w -0.06. (- -), and bT/e 2 w 0.06. ( ), a, c = 0; 

(e)-(f): c = ( ), ce 2 T « 0.035 (- -), and ce 2 T w 0.07. ( ), a, 6 = 0. 



of phase transition without having a significant effect on stability of the dense phase (see 
figures [2aT - [2d|) . The coefficient a has a greater impact on properties of the thin phase 
while b on the properties of the dense phase. At the same time the critical temperature 
increases with c, which affects stability of the dense phase, as shown on figures I2cll2fl 
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Our conclusions on influence of the coefficients a, b and c are based on the properties 
of graphs of a (no) and P(n ) at fixed temperatures. Let us show that they are invariant 
under contact transformations. There are four monomials of generalized degree 5, 
namely, z 7 z 3 , z 3 , z\z± and 2§zs, but only two of them give non-zero Poisson brackets 
with H of generalized degree 10, 

{z 3 z 7 , H } = -3 z 3 zl + 2 z\ (z 7 - z 3 z 4 ) 2 , {z 3 , H } = 2 z% (z 7 - z 3 z 4 ). 

The corresponding 1-parameter groups of contact transformations have the form 

[z 4 ^e Xz3 z 4 , z 5 ^e Xz3 z 5 , z 6 -> e Xza (z 6 + A z 7 ), z 7 -> e A z * z 7 \ , 
[z e ->■ z e + 5 A zl , z 7 -)■ z 7 + A zf] , 

where the omitted variables are left unchanged, and A is a parameter. Both 
transformations preserve the equation of chemical equilibrium fl23|) . since they do 
not change the chemical potentials. The first transformation represents the uniform 
scaling of pressure and concentrations by the scale factor e Xz3 , so it does not affect our 
considerations. The second transformation shifts pressure by the value A z\, which leads 
to violation of equality of P to zero at zero density and hence has no physical sence, so 
the term z| (z 7 — z 3 z 4 ) is restricted in any perturbation of the Debye-Hiickel model. 

6. Conclusion 

As can be seen from the above examples, the introduced method of generalized normal 
forms is a flexible instrument for studying the critical phenomena in various physical 
systems. It can be used whenever there is a natural symplectic or contact form. To 
our knowledge, thermodynamic equations of state have never been considered from 
this point of view. In the context of infinitesimal contact transformations and the 
corresponding systems of ODEs, normal forms were previously studied only in the 
case of quasi-homogeneous unperturbed part H of generalized degree 2 with weight 
7 = (1, . . . , 1, 2) with generic coefficients (see [131 P 37] and references therein). 

In both examples the unperturbed Hamiltonian Hq was a quasi-homogeneous 
polynomial with a canonical weight, and a generalized normal form was obtained rather 
easily. However, it is not so for all interesting systems. For example, the Hamiltonian 



where Ai, A2 > 0, m, n G N, n > m > 1, is not quasi-homogeneous for any canonical 
weight. Such a Hamiltonian arises in studying the interaction between two nonlinear 
oscillators at the level of perturbation. The problem of stability of the equilibrium for a 
system with such an unperturbed part was partially solved in [2] by methods of KAM 
theory: for sufficiently high orders of the perturbation, the origin is Lyapunov stable. At 
low orders it is natural to deal with a normal form instead of an arbitrary perturbation, 
since it contains a significantly smaller number of terms. In the general case, explicit 
formulas for minimal resonant sets turn out to be very complicated. Nevertheless, 



#0 



mp\ + Ai q\ 
2m 



npl + X 2 g| 
2n 
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Lemma 1 gives an algorithm, which can be used to obtain the leading terms of the 
Poisson brackets, and hence a minimal resonant set, for given m, n, N, and linear order 
on monomials. For example, the minimal resonant set for m = 2, n = 3, N = 11, 
and a lexicographical order determined by the order pi >- qi >- p 2 y q 2 consists of 
monomials p 2 q{ x q ] 2 , where 7 < j\ + i 2 + ji < 11 and ji > i 2 , except the following: 

q%qt, qfqh qfqt P2q\q\, P2q\qt, P2q\q\, pIqIqI, p\q\q\, p\q\q\, p\q\q\, plq\qi, 

SS^i 4 4 4^ 4fi 2 4 S 244 24^ S42 S4S S 4 4 

PUlQ^ P2?i5S> P2?i92> P2 9i?2> P2?l92> Pa 9l ?2> P2?l?2> P2?l?2> P2 Ql ?2> P2 ?1 ?2 > 

P% it 12, pi qf qi, p\ qf ql, p\ q\ q\, p\ q\ q 2 , p\ q\ q\, p\ q\ q\, p\ q\ ?2, p\ q\ q\, p\ q\, p\ q\ q%, 

p\ q\ 92, p\ q\ ql, p\ q\, p\ q\ q2, p\ q\, q\ q\, P2 q\ q\, p\ q\ q\, p\ q\ Q2, pi ql, P2 qf qi, p\ q\ q%, 

P% qf, p\ qf - This set consists of 134 monomials, while in an arbitrary perturbation there 
are 1155 terms at these degrees. 
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